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Abstract. Event Shape Data from e + e~ annihilation into hadrons collected by the JADE experiment at 
centre-of-mass energies between 14 GeV and 44 GeV are used to determine the strong coupling as- QCD 
predictions complete to next-to-next-to-leading order (NNLO), alternatively combined with resummed 
next-to-leading-log-approximation (NNLO+NLLA) calculations, are used. The combined value from six 
different event shape observables at the six JADE centre-of-mass energies using the NNLO calculations is 

a s (m z o) = 0.1210 ± 0.0007(stat.) ± 0.0021(exp.) ± 0.0044(had.) ± 0.0036(theo.) 

and with the NNLO+NLLA calculations the combined value is 

a s {m z o) = 0.1172 ± 0.0006(stat.) ± 0.0020(exp.) ± 0.0035(had.) ± 0.0030(theo.) . 

The stability of the NNLO and NNLO+NLLA results with respect to missing higher order contributions, 
studied by variations of the renormalisation scale, is improved compared to previous results obtained with 
NLO+NLLA or with NLO predictions only. The observed energy dependence of as agrees with the QCD 
prediction of asymptotic freedom and excludes absence of running with 99% confidence level. 



1 Introduction 



' Analyses of events originating from e + e annihilation into 
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hadrons allow studies [2 - 5 of Quantum Chromodynamics 
(QCD), the theory of the strong interaction [2H3]. Com- 
parison of observables like jet production rates or event 
shapes with theoretical predictions provides access to the 
determination of the strong coupling as- Recently signifi- 
cant progress in the theoretical calculations of event shape 
observables has been made and next-to-next-to-leading 
order (NNLO) predictions are now available [10] as well 
as matching with resummed calculations in the next-to- 
leading- log-approximation (NLLA) [TT] . As a first applica- 
tion, measurements of as at centre-of-mass-system (cms) 
energies between y/s = 91 GeV and y/s = 206 GeV were 
presented [12]. The same theoretical predictions are used 
in this paper to determine the strong coupling as from 
JADE0 data recorded at lower cms energies. As in [T3] 
and the previous standard LEP and JADE analyses [5] we 
use Monte Carlo simulations to treat hadronisation effects. 
In [13] data for thrust at cms energies 14 < y/s < 207 GeV 
are analysed with combined NNLO+NLLA calculations 
and an analytic model for non-perturbative physics. 



The JADE experiment operated at the PETRA0 e+e" 
collider at DES\q Hamburg, Germany. The data taken 
in the years from 1979 to 1986 cover cms energies between 
12 and 46.4 GeV. 



2 The JADE Detector 

A detailed description of the JADE detector can be found 
in [14]. For this analysis, tracks from charged particles 
and energy deposits in the electromagnetic calorimeter are 
used. The main tracking device, a large volume jet cham- 
ber was located in a 0.48 T solenoidal magnetic field. The 
electromagnetic calorimeter consisted of 2520 lead glass 
blocks in the barrel and 192 lead glass blocks in both end- 
caps with radiation length varying between 9.6 in the end- 
caps and up to 15.7 in the barrel. 



3 Data and Monte Carlo Samples 

The Data and Monte Carlos samples utilised in this anal- 
ysis are identical to those used in a previous determina- 



a The members of the JADE collaboration are listed in [T] 
1 JApan-Deutschland-England 



2 Positron-Elektron- Tandem-Ring- Anlage 

3 Deutsches Elektronen SYnchrotron 
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tion of as [IS] ■ The data correspond to a total integrated 
luminosity of 195/pb taken at average cms energies be- 
tween 14.0 GeV and 43.8 GeV. The breakdown of the 
data samples, including the average cms energies, the en- 
ergy ranges, data taking periods, integrated luminosities 
and the overall numbers of selected hadronic events are 
summarised in table [T] 

Table 1. The average centre-of-mass energy y/s, the energy 
range, data taking period, collected integrated luminosity L 
together with the number of selected data events. 



i/s energy year L selected 

GeV range in GeV (1/pb) events 



14.0 


13.0- 


-15.0 




1981 


1.46 


1783 


22.0 


21.0- 


-23.0 




1981 


2.41 


1403 


34.6 


33.8- 


-36.0 


1981- 


-1982 


61.7 


14313 


35.0 


34.0- 


-36.0 




1986 


92.3 


20876 


38.3 


37.3- 


-39.3 




1985 


8.28 


1585 


43.8 


43.4- 


-46.4 


1984- 


-1985 


28.8 


4374 



Monte Carlo events are generated in large numbers 
to correct the data for experimental acceptance, resolu- 
tion effects and background. Events are simulated using 
PYTHIA 5.7 OH, and for systematic studies with HER- 
WIG 5.9 [17] . Subsequently the generated events are pro- 
cessed through a full simulation of the JADE detector and 
are reconstructed in the same way and with the same pro- 
gram chain as the data. For comparison with the corrected 
data and for the correction of hadronisation effects large 
samples of Monte Carlo events have been produced using 
PYTHIA 6.158, HERWIG 6.2 and ARIADNE 4.11 [TS]. 
We use the model parameters as determined at -^=91 GeV 
by the OPAL experiment at the LEP e+e" collider [T9T2T)] . 



4 Experimental Procedure 

4.1 Event Selection 

The selection of identified and well measured hadronic 
event candidates follows the procedure outlined in [TS] . 
Events with a large momentum imbalance due to photons 
emitted in the initial state are rejected. The event selec- 
tion is based on minimal requirements for charged particle 
multiplicity, visible energy and longitudinal momentum 
imbalance. The dominating backgrounds from hadronic r 
decays and two-photon interactions with hadronic final 
states are supressed to negligble levels. 

4.2 Event Shape Distributions 

The properties of hadronic events can be described by 
event shape observables. Event shape observables used 
for this analysis are thrust (1 — T) [2TJ[22], heavy jet 
mass (Mh) [23], wide and total jet broadening (i?w and 
B T ) [23], C-Parameter (C) [25H27] and the transition value 



between 2 and 3 jet configurations [28][25] defined by the 
Durham jet algorithm (y^ 3 ) [50] . Whenever we refer to a 
generic event shape observable 1 — T, Mh, Bt, B^, C or 
IJ23 we use the symbol y. 

The event shape observables are calculated from se- 
lected charged particle tracks and calorimeter clusters af- 
ter correcting for double counting of energy as described 
in [TS]. We compared the data with the predictions from 
Monte Carlo simulations as described above and found 
good agreement at all energy points. Similar observations 
were made in [TS] for related observables. 

4.3 Corrections to the Data 

The corrections to the data for limited experimental res- 
olution, acceptance and bb background follows exactly 
the treatment described in [15]. Selected charged parti- 
cle tracks as well as electromagnetic clusters are used to 
calculate the event shape observables. 

From simulated events two different distributions are 
built: the detector-level distribution and the hadron-level 
distribution. The detector-level distributions are calculated 
exactly in the same way as for data using measured charged 
particle tracks and calorimeter clusters. The hadron-level 
distributions use the true four-momenta of the stable par- 
ticle^ in events where the centre-of-mass energy is re- 
duced due to initial state radiation (ISR) by less than 
0.15 GeV. The bin- by-bin corrections for the data dis- 
tributions are derived from the ratio of hadron-level to 
detector-level distributions for simulated events with u, 
d, s or c primary quarks. Contributions from B hadron 
decays bias the measurement of event shape observables 
and therefore the expected contribution from e + e~ — > bb 
events is subtracted from the detector-level distributions 
before the corrections are applied. The simulations were 
optimised by OPAL to describe production and decays 
of B hadrons [HJ[2n] . The good description of our uncor- 
rected data by the simulations confirms that using the 
simulations to subtract the e + e _ —> bb background is jus- 
tified. 

In order to study systematic uncertainties the selec- 
tion and correction procedures are modified and the whole 
analysis is repeated. The evaluation of the systematic un- 
certainties follows identically the procedure described in [TS] 



5 Determinations of as 
5.1 QCD Calculations 

The distributions of the event shape observables are pre- 
dicted by (NNLO) perturbative QCD calculations [10]: 

1 da dA dB 9 dC , 

-~r = -r a s + -rcii + —a-Z (1) 

<7dy 0(a 3 s) dy dy dy 

4 All particles with lifetimes greater than 300 ps are consid- 
ered stable. 
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with as = as(/i)/(27r). Equation ([T]) is shown for renor- 
malisation scale /i — Q, where Q is the physical scale usu- 
ally identified with the cms energy y/s for hadron produc- 
tion in e + e~ annihilation. The coefficient distributions for 
leading order (LO) dA/dy, next-to-leading order (NLO) 
dB/dy and NNLO dC/dy were kindly provided by the 
authors of [10]. In [31] a problem at small values of y 
with the NNLO terms calculated in [TU] was shown, but 
it does not affect the kinematic regions selected in our 
fits. The normalisation to the total hadronic cross sec- 
tion and the terms generated by variation of the renor- 
malisation scale parameter = /i/Q are implemented 
according to [IDj . The prediction in equation JT]) may be 
combined with resummed NLLA calculations [TT| using 
the In i?-matching scheme; we refer to these predictions as 
NNLO+NLLA. The In i?-matching procedure ensures that 
in the combination of the fixed order NNLO calculations 
and the resummed NLLA calculations no double counting 
of common terms occurs. The NNLO+NLLA predictions 
are here compared with experimental data for the first 
time. 

The theoretical predictions provide distributions at the 
level of quarks and gluons, the so-called parton-level. The 
distributions calculated using the final state partons af- 
ter termination of the parton showering in the models 
are also said to be at the parton-level. To compare the 
QCD predictions with measured hadron-level event shape 
distributions the predictions are corrected for hadroni- 
sation effects. These corrections are obtained by calcu- 
lating in the Monte Carlo models the ratio of the cu- 
mulative distributions at hadron-level and parton-level. 
The corrections are applied to the cumulative prediction 
R(y) = J y 1/ada/dy'dy' as in [32]. 

It was shown in [T2] that the event shape observable 



ered by the combined hadronisation and theory systematic 
variations. 



5.2 Measurements of as 

The strong coupling as is determined by a \ 2 fit to each 
of the measured event shape distributions at the hadron- 
level, i.e. corrected for experimental effects. A x 2 value is 
calculated at each cms energy: 



(d t -U(a s ))(V-%(d ] -t 3 (a s )) (2) 



where i,j count the bins within the fit range of the event 
shape distribution, di is the measured value in the iih bin, 
^i(as) is the QCD prediction for the ith bin corrected for 
hadronisation effects, and Vij is the covariance matrix of 
the di. The final prediction is obtained by integrating the 
QCD predictions in equation ([1]) over the bin width after 
application of the hadronisation correction as explained 
above in section 15.11 The x 2 value is minimised with re- 
spect to as while the renormalisation scale factor is set to 
Xfj, = 1. 

The evolution of the strong coupling as as a function 
of the renormalisation scale is implemented in three loops 
as shown in [5]. Since the evolution of as in the cms en- 
ergy range considered here does not involve the crossing 
of flavour thresholds it does not introduce significant un- 
certainties. In order to quantify the uncertainty from the 
evolution procedure we evolve as(mzo) = 0.118 from myp 
to 14 GeV in three loops and two loops and find a relative 
difference of 0.1%. 

In order to take the correlations between different bins 



distributions derived from the parton-level of the Monte into account the covariance matrix V l3 is computed fol- 
Carlo generators are described reasonably well by the NNLO lowin S the approach described in [53]: 



calculation in their fit ranges. We compared the parton- 
level predictions of the Monte Carlo generators with the 
QCD predictions in NNLO or NNLO+NLLA with a s (m z o ) = 

0. 118 and x M = 0.5, 1.0, 2.0 at y/s = 14, 22, 35 and 44 GeV. 
We study the quantity r(y) theoryiM o = da / dy theory /da / dy M c 

1. In addition we compute the corresponding quantities 
r (y)MCi.MCj for any pair i,j of Monte Carlo predictions 
and r(y)a: M= o.5;i;2,x M =o.5;i;2 for theory predictions with dif- 
ferent renormalisation scale values. The largest values of 

the ahs(r(y) M C i ,MC j ) and the abs(r(y) :EM= o.5;i;2,x M =o.5;i;2) 
at each y are added in quadrature to define the uncer- 
tainty Ar(y) of r(y)theory,MC ■ The average values r of the 
abs(r(y)) over the fit ranges (see below) are taken as a 
measure of the consistency between between theory and 
Monte Carlo predictions. The ratios of f theory, MC with the 
average error Ar(y) are generally smaller than or about 
equal to unity and reach values of about two for C at 
y/s = 44 GeV. The model dependence of the hadronisa- 
tion correction and the renormalisation scale dependence 
of the theory will be studied as systematic uncertainties 
below. Our studies show that systematic uncertainties in- 
troduced by discrepancies between the theory predictions 
and the Monte Carlo parton-level predictions will be cov- 



Vn = 



dPj dPj 
dN k dN k 



E 

k 



N k 



(3) 



N5 ik - N, 



i) (NS jk 



Ni is the number of data events, iVj = cti(Ni — &,) the 
number of events after subtraction of background bi from 
bb events and multiplying by a correction on for detector 
effects, Pi is the normalised hadron-level distribution at 
bin i and N = J2 k Nk- 

The fit ranges are determined by several considera- 
tions. We require the leading log terms to be less than 
unity, because we also use a fixed order expansion without 
rcsummation of log-enhanced terms, see e.g. [34) . and for 
consistency we use the same fit ranges in the NNLO and 
NNLO+NLLA analyses. The leading log term of dA/dy is 
hiy/y and we require dslny/z/ < 1 for as(m z o) = 0.118, 
14 < \fs < 200 GeV and y > 0.1. The upper limit is given 
by the requirement that all three orders of the NNLO cal- 
culations contribute, i.e. the fit range extends to the kine- 
matic limit of the LO coefficients dA/dy. The fit range for 
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Mh should be compared with the related fit range for 1— T 
after squaring, because My ~ 1 — T in LO. The fit ranges 
for C and 1 — T are related by a factor of (ln6)/6 ~ 0.3. 
The resulting fit ranges are shown in table [5] 



Table 2. Fit ranges at all cms energies 



1 — T 
0.10-0.27 



0.26-0.50 



St 
0.16-0.30 



Bw 
0.10-0.23 



C 

0.34-0.72 



D 
2/23 

0.01-0.20 



The detector corrections (see section |4~3]) are generally 
±20% or less within the fit ranges. The hadronisation cor- 
rection factors are maximally 3.5 for Bt at yfs = 14 GeV 
but are generally smaller than 2 for the larger cms en- 
ergy points. The hadronisation corrections are smallest 
and have the least variations over the fit ranges for yJ? 3 . 

The evaluation of the systematic errors of the as mea- 
surements takes into account experimental effects, the hadro-inate the total errors (except for M H and C), because at 
nisation correction procedure and uncertainties of the the- 



5.3 Results from NNLO Fits 



The results of the NNLO fits are summarised in table [3] 
In figure Q] the 1 — T event shape distributions together 
with the NNLO fit results for the six energy points are 
shown. The % 2 /d.o.f. values are between 0.7 for yfs = 
14 GeV and 2.5 for yfs = 34.6 GeV. The fit results for the 
other event shape observables return similar results with 
0.3 < x 2 /d.o.f. < 3.8. We note that at cms energies where 
we have big data samples the % 2 /d-o.f. value tend to be 
larger. The % 2 /d.o.f. values are based on statistical errors 
only while the combined experimental and hadronisation 
uncertainties are at least a factor of two larger than the 
statistical errors leading to a reduction of x 2 /d.o.f. by a 
factor of at least 4 if these uncertainties were taken into 
account in the fits. We conclude that there is no significant 
disagreement between the event shape data and the QCD 
fits. 

The results at each cms energy are remarkably con- 
sistent with each other, we find root-mean-square (rms) 
values for a s (y/s) between 0.003 at 44 GeV and 0.008 at 
22 GeV. The hadronisation uncertainties at 14 GeV dom- 



14 GeV the hadronisation corrections are largest. The sta- 
tistical errors are sizeable at yfs = 14, 22, 38 and 44 GeV, 
where there is only limited statistics and quite small at 
yfs = 34.6 and 35 GeV where we have large data samples. 
The experimental uncertainties depend somewhat on the 
cms energy with smaller values at higher yfs where we 
have larger data samples. 

The hadronisation uncertainties for C are the largest 
except at yfs = 14 GeV. This as been observed before [3"6l - 



ory. The three sources of systematic uncertainty are added 
in quadrature to the statistical error taken from the fits 
to obtain the total errors. Below we describe how we find 
the systematic uncertainties: 

Experimental Uncertainties The analysis is repeated with 
slightly varied event and track selection cuts and a sys- 
tematic uncertainty from variation of the fit ranges is 
studied [TS]. The cross section used in the subtraction 
of bb events is varied by ±5% which takes account 
of possible differences in the efficiency determination 
using the simulations of e + c~ — > qq (q=u,d,s,c) and 
e + e _ — > bb events. For each experimental variation 
the value of as is determined and compared to the 
central (default) value. The quadratic sum of the dif- 
ferences and the fitrange uncertainty is taken as the 
experimental systematic uncertainty. 

Hadronisation For the default analysis, PYTHIA is used 
to estimate the corrections originating from hadroni- 
sation effects (section l5.1[) . As a systematic variation 
HERWIG and ARIADNE are used to evaluate the ef- 
fects of hadronisation. The larger of the deviations is 
taken as systematic hadronisation uncertainty. It was 
observed in [35K2S] that systematic uncertainties be- 
tween the PYTHIA, HERWIG and ARIADNE models 
are generally much larger than systematic uncertain- 
ties from varying the parameters of a given model. 

Theoretical Uncertainties The theoretical prediction of eventu nce rtainties of the NNLO+NLLA~ fits are also smaller 



5.4 Results from NNLO+NLLA Fits 



The results of the NNLO+NLLA fits are given in tabled] 
The values of % 2 /d.o.f. are slightly smaller than the corre- 
sponding values of the NNLO fits in most cases, indicating 
a somewhat better consistency of the NNLO+NLLA cal- 
culations with the fitted data points. The rms values of 
a s (yfs) are between 0.007 at yfs = 22 GeV and 0.003 at 
yfs = 44 GeV, i.e. the scatter of individual results is es- 
sentially the same as for the NNLO analysis. The pattern 
of statistical errors and experimental and hadronisation 
uncertainties is the same as for the NNLO fits discussed 
above. Compared with the NNLO analysis the theoreti- 
cal uncertainties are reduced by 10 — 20% and the val- 
ues of as are lower by 4% on average. The hadronisation 



shape observables is a finite power series in as • The un- 
certainties originating from missing higher order terms 
are assessed by changing the renormalisation scale fac- 
tor to x^ — 0.5 and = 2.0. The larger deviation 
from the default value of as is taken as the systematic 
uncertainty. 



in most cases. The difference in as between NNLO and 
NNLO+NLLA calculations is smaller than the difference 
in as between NLO and NLO+NLLA calculations as ex- 
pected in [11]. As discussed above in section l5~3l for NNLO 
fits the sometimes large x 2 /d.o.f. values can be explained 
by the small statistical errors in some data sets. 
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Table 3. Results of NNLO fits to event shape observable distributions at the JADE cms energies. 



[GeV] 


Obs. 


as(Vs) 


±stat. 


±exp. 


±had. 


±theo. 


X 2 /d.o.f. 


f4.0 


1 — T 


0.1587 


0.0098 


0.0213 


0.0366 


0.0100 


4.0/6 


14.0 


M H 


0.1759 


0.0080 


0.0133 


0.0093 


0.0099 


9.0/5 


14. u 


D 

n>T 


U.lOo/ 


U.UUoO 


U.uuyo 


n OQQ7 


U.Uloz 


1.0/0 


14.0 


Bw 


0.1730 


0.0053 


0.0088 


0.0188 


0.0088 


3.8/5 


14.0 


C 


0.1583 


0.0150 


0.0169 


0.0113 


0.0089 


9.2/5 


14.0 


y? 3 


0.1671 


0.0039 


0.0054 


0.0101 


0.0063 


11.5/8 


22.0 


1 - T 


0.1410 


0.0075 


0.0054 


0.0195 


0.0070 


6.8/6 


22.0 


Mn 


0.1555 


0.0070 


0.0105 


0.0064 


0.0061 


5.3/5 


22.0 


Bt 


0.1399 


0.0070 


0.0050 


0.0178 


0.0076 


4.0/5 


22.0 


Bw 


0.1551 


0.0046 


0.0046 


0.0106 


0.0060 


10.2/5 


22.0 


C 


0.1385 


0.0084 


0.0073 


0.0332 


0.0062 


3.7/5 


22.0 


!& 


0.1545 


0.0033 


0.0025 


0.0084 


0.0049 


10.6/8 



34.6 


1 - T 


0.1396 


0.0017 


0.0040 


0.0100 


0.0069 


14.8/6 


34.6 


M H 


0.1477 


0.0017 


0.0070 


0.0025 


0.0053 


18.9/5 


34.6 


Bt 


0.1392 


0.0016 


0.0025 


0.0063 


0.0076 


6.4/5 


34.6 


Bw 


0.1457 


0.0013 


0.0047 


0.0049 


0.0049 


3.5/5 


34.6 


C 


0.1374 


0.0017 


0.0041 


0.0126 


0.0062 


12.7/5 


34.6 


2/23 


0.1404 


0.0009 


0.0012 


0.0061 


0.0035 


14.3/8 


35.0 


1 — T 


0.1430 


0.0014 


0.0040 


0.0094 


0.0074 


14.0/6 


35.0 


Mn 


0.1532 


0.0014 


0.0065 


0.0023 


0.0059 


7.0/5 


35.0 


Bt 


0.1432 


0.0013 


0.0042 


0.0064 


0.0082 


13.8/5 


35.0 


Bw 


0.1496 


0.0011 


0.0063 


0.0046 


0.0054 


10.7/5 


35.0 


C 


0.1427 


0.0014 


0.0036 


0.0118 


0.0069 


9.3/5 


35.0 


D 
2/23 


0.1451 


0.0008 


0.0021 


0.0061 


0.0039 


15.0/8 



38.3 


1 — T 


0.1427 


0.0049 


0.0064 


0.0081 


0.0073 


11.3/6 


38.3 


M H 


0.1578 


0.0048 


0.0067 


0.0022 


0.0066 


2.0/5 


38.3 


Bt 


0.1447 


0.0045 


0.0068 


0.0048 


0.0085 


1.8/5 


38.3 


Bw 


0.1488 


0.0039 


0.0069 


0.0034 


0.0053 


3.1/5 


38.3 


C 


0.1369 


0.0049 


0.0036 


0.0101 


0.0061 


9.7/5 


38.3 


D 
2/23 


0.1380 


0.0028 


0.0038 


0.0071 


0.0032 


23.5/8 


43.8 


1 — T 


0.1341 


0.0029 


0.0034 


0.0057 


0.0061 


11.3/6 


43.8 


Mn 


0.1403 


0.0029 


0.0062 


0.0014 


0.0046 


11.5/5 


43.8 


Bt 


0.1312 


0.0027 


0.0027 


0.0043 


0.0064 


9.8/5 


43.8 


Bw 


0.1373 


0.0023 


0.0050 


0.0028 


0.0041 


10.7/5 


43.8 


C 


0.1342 


0.0028 


0.0042 


0.0083 


0.0058 


4.3/5 


43.8 


2/2°3 


0.1333 


0.0017 


0.0026 


0.0054 


0.0029 


21.0/8 



5.5 Combination of Results 



The results obtained at each energy point for the six event 
shape observables are combined using error weighted av- 
eraging as in [5 , 33 , 39 , 40 . The statistical correlations be- 
tween the six event shape observables are estimated at 
each energy point from fits to hadron-level distributions 
derived from 50 statistically independent Monte Carlo sam- 
ples. The experimental uncertainties are determined as- 
suming that the smaller of a pair of correlated experimen- 
tal errors gives the size of the fully correlated error (partial 
correlation) . The hadronisation and theory systematic un- 
certainties are found by repeating the combination with 
changed input values, i.e. using a different hadronisation 
model or a different value of x^. The results are given in 
table [5] and shown for the NNLO analysis in figure [51 be- 



cause this allows a direct comparison with the results of 
the NNLO analysis of ALEPH event shape data . 

The statistical uncertainties of the combined results 
are reduced as expected. The systematic uncertainties of 
the combined results tend to be close to the best values 
from individual observables, because the systematic un- 
certainties are not completely correlated and because ob- 
servables with smaller uncertainties have larger weights in 
the combination procedure. 

A combination of the combined results at the six JADE 
energy points shown in table [5] after running to a common 
reference scale m^o using the combination procedure de- 
scribed above results in 

a s {m z o) = 0.1210 ±0.0007(stat.) ± 0.0021(cxp.) (4) 
±0.0044(had.) ± 0.0036(theo.) 
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Table 4. Results of NNLO+NLLA fits to event shape observable distributions at the JADE cms energies. 



[GeV] 


Obs. 


as(Vs) 


±stat. 


±exp. 


±had. 


±theo. 


X 2 /d.o.f. 


14.0 


1 — T 


0.1543 


0.0092 


0.0197 


0.0362 


0.0089 


4.3/6 


14.0 


M H 


0.1641 


0.0065 


0.0105 


0.0126 


0.0095 


9.2/5 


14.0 


Bt 


0.1620 


0.0078 


0.0083 


0.0343 


0.0122 


1.5/5 


14.0 


Bw 


0.1540 


0.0038 


0.0060 


0.0157 


0.0072 


1.4/5 


14.0 


C 


0.1466 


0.0128 


0.0151 


0.0131 


0.0063 


11.1/5 


14.0 


y? 3 


0.1661 


0.0038 


0.0057 


0.0129 


0.0060 


8.0/8 


22.0 


1 - T 


0.1383 


0.0070 


0.0051 


0.0185 


0.0061 


6.1/6 


22.0 


Mn 


0.1464 


0.0058 


0.0089 


0.0040 


0.0049 


4.8/5 


22.0 


Bt 


0.1350 


0.0064 


0.0044 


0.0159 


0.0062 


3.4/5 


22.0 


Bw 


0.1408 


0.0034 


0.0035 


0.0074 


0.0047 


7.2/5 


22.0 


C 


0.1337 


0.0078 


0.0064 


0.0323 


0.0061 


4.3/5 


22.0 


y% 


0.1538 


0.0033 


0.0022 


0.0090 


0.0045 


8.5/8 



34.6 


1 - T 


0.1365 


0.0016 


0.0036 


0.0092 


0.0057 


13.7/6 


34.6 


M H 


0.1399 


0.0014 


0.0055 


0.0019 


0.0041 


15.9/5 


34.6 


Bt 


0.1338 


0.0014 


0.0021 


0.0055 


0.0058 


6.3/5 


34.6 


Bw 


0.1332 


0.0010 


0.0034 


0.0037 


0.0035 


6.2/5 


34.6 


C 


0.1326 


0.0016 


0.0037 


0.0111 


0.0061 


6.1/5 


34.6 


2/23 


0.1401 


0.0009 


0.0011 


0.0059 


0.0030 


5.5/8 



35.0 


1 — T 


0.1399 


0.0013 


0.0037 


0.0086 


0.0062 


10.5/6 


35.0 


Mn 


0.1444 


0.0012 


0.0050 


0.0015 


0.0047 


8.1/5 


35.0 


Bt 


0.1376 


0.0012 


0.0038 


0.0055 


0.0065 


11.9/5 


35.0 


Bw 


0.1363 


0.0008 


0.0048 


0.0035 


0.0040 


12.8/5 


35.0 


C 


0.1373 


0.0013 


0.0032 


0.0104 


0.0068 


4.6/5 


35.0 


D 
2/23 


0.1447 


0.0008 


0.0022 


0.0059 


0.0035 


9.6/8 



38.3 


1 — T 


0.1400 


0.0046 


0.0060 


0.0073 


0.0064 


10.2/6 


38.3 


M H 


0.1484 


0.0040 


0.0050 


0.0015 


0.0056 


1.9/5 


38.3 


Bt 


0.1390 


0.0041 


0.0060 


0.0043 


0.0068 


1.3/5 


38.3 


Bw 


0.1360 


0.0030 


0.0052 


0.0026 


0.0042 


2.5/5 


38.3 


C 


0.1329 


0.0045 


0.0034 


0.0085 


0.0064 


7.4/5 


38.3 


D 
2/23 


0.1385 


0.0028 


0.0040 


0.0063 


0.0034 


18.1/8 


43.8 


1 — T 


0.1313 


0.0027 


0.0031 


0.0052 


0.0048 


11.6/6 


43.8 


Mn 


0.1337 


0.0025 


0.0049 


0.0011 


0.0034 


10.8/5 


43.8 


Bt 


0.1265 


0.0025 


0.0023 


0.0039 


0.0047 


10.0/5 


43.8 


Bw 


0.1269 


0.0018 


0.0037 


0.0022 


0.0027 


11.1/5 


43.8 


C 


0.1295 


0.0025 


0.0037 


0.0072 


0.0054 


3.5/5 


43.8 


2/2°3 


0.1331 


0.0017 


0.0024 


0.0050 


0.0024 


18.6/8 



(a s (m z o) = 0.1210 ± 0.0061) for the NNLO analysis and 

a s (m z o) = 0.1172 ±0.0006(stat.) ± 0.0020(cxp.) (5) 
±0.0035(had.) ± 0.0030(theo.) 

(a s (m z o) = 0.1172 ± 0.0051) for the NNLO+NLLA anal- 
ysis. The NNLO+NLLA result has smaller hadronisation 
and theory uncertainties compared with the values in the 



for ag( m z°) °f 4% is among the most precise determina- 
tions of as(?7iz ) currently available. 

After running the fit results for as{y/s) 101 each ob- 
servable to the common reference scale m z o we combine 
the results for a given observable to a single value. We 
use the same method as above and obtain the results for 
as(m z o) shown in table [5] The rms values of the results 
NNLA analysis. We choose the latter result from NNLO+NLIfoV a s (m z o) are 0.0029 for the NNLO analysis and 0.0026 
fits as our final result, because it is based on the most com- for the NNLO+NLLA analysis; both values are consis- 
plete theory predictions and it has smaller theory uncer- tent with the errors of the corresponding combined results 
tainties. It is consistent with the world average of as (m z o) = shown in equations (|U) and (JSJ) - Figure [3] shows the com 
0.119+0.001 [H], the recent NNLO analysis of event shape 
data from the ALEPH experiment as(m z o) — 0.1240 + 
0.0033 [12] as well as with the related average of as (wz° ) = 
0.120 ± 0.005 from the analyses of the LEP experiments 
using NLO+NLLA QCD predictions [5]. The total error 



bined results of as(m z o) for each observable together with 
results from alternative analyses discussed below. Com- 
bining the combined results for each observable or com- 
bining all individual results after evolution to the com- 
mon scale m z o yields results consistent with equation ([SJ 
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Fig. 1. The plots show as points with statistical error bars the 
1 — T distributions at hadron level at y/s = 14 to 43.8 GeV. 
Some error bars are smaller than the data points. Superim- 
posed as histograms are the NNLO+NLLA predictions com- 
bined with hadronisation effects using the corresponding fit 
results for as(y/s) shown in table [3] The arrows indicate the 
fit ranges. 

Table 5. Combined values of as(i/s) at the JADE cms en- 
ergies from NNLO (upper section) and NNLO+NLLA (lower 
section) analyses together with the statistical, experimental, 
hadronisation and theory errors. 



>/S [GeV] 




±stat. 


±exp. 


±had. 


±theo. 


14.0 


0.1690 


0.0046 


0.0065 


0.0124 


0.0076 


22.0 


0.1527 


0.0040 


0.0036 


0.0090 


0.0056 


34.6 


0.1420 


0.0012 


0.0025 


0.0058 


0.0050 


35.0 


0.1463 


0.0010 


0.0032 


0.0059 


0.0055 


38.3 


0.1428 


0.0033 


0.0045 


0.0060 


0.0051 


43.8 


0.1345 


0.0021 


0.0031 


0.0043 


0.0045 


14.0 


0.1605 


0.0044 


0.0065 


0.0148 


0.0073 


22.0 


0.1456 


0.0036 


0.0033 


0.0077 


0.0048 


34.6 


0.1367 


0.0011 


0.0023 


0.0046 


0.0040 


35.0 


0.1412 


0.0009 


0.0032 


0.0049 


0.0047 


38.3 


0.1388 


0.0030 


0.0043 


0.0042 


0.0048 


43.8 


0.1297 


0.0019 


0.0028 


0.0033 


0.0034 



0.2 

S 
0.19 

0.18 

0.17 

0.16 

0.15 

0.14 

0.13 

0.12 

0.11 

0.1 



JADE NNLO 

■ ALEPH NNLO 
1-T, Mjj, By, 15 yy, C, 




a s (m 7 )=0.1210±0.0061 



Vs [GeV] 



Fig. 2. The values for as at the JADE energy points. The 
inner error bars correspond to the combined statistical and 
experimental errors and the outer error bars show the total 
errors. The results from y^s = 34.6 and 35 GeV have been 
combined for clarity. The full and dashed lines indicate the 
result from our JADE NNLO analysis as shown on the figure. 
The results from the NNLO analysis of ALEPH data [12] are 
shown as well. 

Table 6. Combined values of as (m z o ) for each observable from 
NNLO (upper section) and NNLO+NLLA (lower section) anal- 
yses together with the statistical, experimental, hadronisation 
and theory errors. 



Obs. 


cts(m z o) 


±stat. 


±exp. 


±had. 


±theo. 


1-T 


0.1196 


0.0011 


0.0028 


0.0067 


0.0049 


M H 


0.1266 


0.0009 


0.0047 


0.0014 


0.0040 


Bt 


0.1190 


0.0009 


0.0023 


0.0047 


0.0055 


Bw 


0.1232 


0.0008 


0.0034 


0.0037 


0.0035 


C 


0.1184 


0.0013 


0.0029 


0.0081 


0.0045 


!& 


0.1201 


0.0005 


0.0014 


0.0046 


0.0026 


1-T 


0.1175 


0.0010 


0.0026 


0.0061 


0.0041 


M H 


0.1210 


0.0008 


0.0037 


0.0011 


0.0032 


Bt 


0.1151 


0.0009 


0.0019 


0.0039 


0.0042 


Bw 


0.1143 


0.0006 


0.0026 


0.0028 


0.0026 


C 


0.1148 


0.0011 


0.0027 


0.0073 


0.0044 


J/2 D 3 


0.1199 


0.0005 


0.0013 


0.0046 


0.0023 



within Aas{m z a) = 0.0004 and the uncertainties also 
agree. 

The hadronisation uncertainty of Mh at each energy 
point and in the combinations shown in table [5] is the 
smallest. We have repeated the combinations without Mh 
and found results for asijnzo) consistent within 0.6% with 
our main results with hadronisation uncertainties increased 
by 14% (NNLO) or 20% (NNLO+NLLA). 



In order to study the compatibility of our data with the 
QCD prediction for the evolution of the strong coupling 
with cms energy we repeat the combinations with or with- 
out evolution of the combined results to the common scale. 
We set the theory uncertainties to zero since these uncer- 
tainties are highly correlated between energy points. We 
conservatively assume the hadronisation uncertainties to 
be partially correlated, because these uncertainties depend 
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NNLO 



NNLO 
+NLLA 



NLO 
+NLLA 



NLO 



o.i 



Fig. 3. The combined results for as(m z o) for each type of anal- 
ysis as indicated on the figure. The shaded bands and dashed 
lines show the combined values of as(m z o) with total errors. 
The inner and outer error bars show the combined statistical 
and experimental and the total errors. 



strongly on the cms energy. The x 2 probabilities of the av- 
erages for running (not running) with NNLO+NLLA fits 
then become 0.39 (9.9 • 10~ 3 ). With the NNLO fits the X 2 
probabilities for running (not running) are 0.48 (1.2-10~ 3 ). 
We interpret this as strong evidence for the dependence of 
the strong coupling on cms energy as predicted by QCD 
from JADE data alone. 



larger by about 60%. The analysis using NLO predic- 
tions gives theoretical uncertainties larger by a factor of 
2.6 and the value for as(wzo) is larger compared to the 
NNLO or NNLO+NLLA results. It has been observed 
previously that values for as from NLO analysis with 
Xfj, = 1 are large in comparison with most other analy- 
ses [44] . Both the NLO+NLLA or NNLO+NLLA anal- 
yses yield a smaller value of as{m z o) compared to the 
respective NLO or NNLO results. The difference between 
NNLO+NLLA and NNLO is smaller than the difference 
between NLO+NLLA and NLO, since a larger part of the 
NLLA terms is included in the NNLO predictions. 



5.7 Renormalisation Scale Dependence 

The theoretical uncertainty due to missing higher order 
terms is evaluated by setting the renormalisation scale pa- 
rameter x^ to 0.5 or 2.0. In order to assess the dependence 
of as on the renormalisation scale the fits are repeated 
using NNLO, NNLO+NLLA, NLO and NLO+NLLA pre- 
dictions with 0.1 < < 10. The strong coupling as 
as well as the x 2 /d.o.f. as a function of x^ for 1 — T 
at y/s = 35 GeV are shown in figure |H The x 2 /d.o.f. 
curves for the NLO+NLLA and NNLO+NLLA fits show 
no local minimum in the x^ range studied. The as(mzo) 
values from NLO predictions are the largest for x M > 0.2. 
The a s (m z o) values using NLO+NLLA and NNLO cal- 
culation almost cross at the natural choice of the renor- 
malisation scale x^ = 1 while the as(m.z°) value from 
the NNLO+NLLA fit is slightly lower. The NLLA terms 
are almost identical to the 0(a|) -terms in 
the NNLO calculation. A similar behaviour can be ob- 
served for Bt and y^- The slopes of the as(mzo) curves 
of the NNLO and NNLO+NLLA fits around the default 
choice X/j, = 1 are smaller than the slopes for the NLO 
and NLO+NLLA fits leading to the decreased theoretical 
uncertainties in our analyses. 



5.6 Comparison with NLO and NLO+NLLA 

For a comparison of our results with previous as mea- 
surements the fits to the event shape distributions are re- 
peated with NLO predictions and with NLO predictions 
combined with resummed NLLA with the modified In R- 
matching scheme (NLO+NLLA), both with x^ = 1. The 
NLO+NLLA predictions with the modified In i?-matching 
scheme were the standard of the final analysis of the LEP 
experiments [33,40,42,43 . The fit ranges as well as the 
procedures for evaluation of the systematic uncertainties 
are identical to the ones in our NNLO and NNLO+NLLA 
analyses. 

The combination of the fits using NLO predictions re- 
turns a s (m z o) = 0.1301 ± 0.0009(stat.) ± 0.0029(exp.) ± 
0.0054(had.)±0.0086(theo.), the fits using combined NLO+ 
NLLA predictions yield a s (m z o) = 0.1172±0.0007(stat.)± 
0.0022(exp.) ± 0.0039(had.) ± 0.0054(theo.) and these re- 
sults are shown in figure [3l The result obtained with the 
NLO+NLLA prediction is consistent with our NNLO and 
NNLO+NLLA analyses, but the theory uncertainties are 



6 Conclusion 

In this paper we present measurements of the strong cou- 
pling as using event shape observable distributions at cms 
energies 14 < y/s < 44 GeV. To determine as fits us- 
ing NNLO and combined NNLO+NLLA predictions were 
used. Combining the results from NNLO+NLLA fits to 
the six event shape observables 1 — T, Mh, i?W: Bt, C 
and yp 3 at the six JADE energy points returns as(?nz°) = 
0.1172±0.0006(stat.)±0.0020(exp.)±0.0035(had.)±0.0030(theo.), 
with a total error on as(m z o) of 4%. The investigation of 
the renormalisation scale dependence of as (mzo ) shows a 
reduced dependence on x^ when NNLO or NNLO+NLLA 
predictions are used, compared to analyses with NLO or 
NLO+NLLA predictions. The more complete NNLO or 
NNLO+NLLA QCD predictions thus lead to smaller the- 
oretical uncertainties in our analysis. The combined re- 
sults for «s at each cms energy are consistent with the 
running of as as predicted by QCD and exclude absence 
of running with a confidence level of 99%. 



S. Bethke et al.: Determination of the Strong Coupling as 



9 




0.11 - 

J _l I I I I I 

4 _ NNLO+NLLA 

\ NNLO 

\ NLO+NLLA 

3 - \ NLO 



2 




J U i i I i i I 

10 1 1 10 



Fig. 4. The plot shows the result of as(m z o) and x 2 /d.o.l. of 
the fit to the thrust event shape distribution for \/s = 35 GeV. 
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